molmass_acetate = 98; % Molar mass of potassium acetate
close all;
%% 37C acetate
tempmat = monod_maxOD_acetate_37C;
smat = size(tempmat);

conc1 = tempmat([1 2 3 4],1);
odmax1 = OD_calibrated_blanked(nanmean(tempmat([1 2 3 4],2:4),2));
odmax_err1 = OD_calibrated_blanked(nanstd(tempmat([1 2 3 4],2:4),0,2));

% Linear scale
figure;
errorbar(conc1/100*1000/molmass_acetate*1000,odmax1,odmax_err1, 'ro-', 'linewidth', 1);

set(gcf, 'position', [0 0 400 300])
set(gca, 'fontsize', 20);
xlabel('Concentration (mM)')
ylabel('Max OD600')
box off;
xlim([0 125]);


%% 37C acetate
tempmat = monod_maxOD_acetate_30C;
smat = size(tempmat);

conc1 = tempmat([1 2 3 4],1);
odmax1 = OD_calibrated_blanked(nanmean(tempmat([1 2 3 4],2:4),2));
odmax_err1 = OD_calibrated_blanked(nanstd(tempmat([1 2 3 4],2:4),0,2));

% Linear scale
hold on;
errorbar(conc1/100*1000/molmass_acetate*1000,odmax1,odmax_err1, 'ko-', 'linewidth', 1);

set(gcf, 'position', [0 0 400 300])
set(gca, 'fontsize', 20);
xlabel('Concentration (mM)')
ylabel('Max OD600')
box off;
xlim([0 125]);


%% 37C acetate
tempmat = monod_maxOD_acetate_25C;
smat = size(tempmat);

conc1 = tempmat([1 2 3 4],1);
odmax1 = OD_calibrated_blanked(nanmean(tempmat([1 2 3 4],2:4),2));
odmax_err1 = OD_calibrated_blanked(nanstd(tempmat([1 2 3 4],2:4),0,2));

% Linear scale
hold on;
errorbar(conc1/100*1000/molmass_acetate*1000,odmax1,odmax_err1, 'bo-', 'linewidth', 1);

set(gcf, 'position', [0 0 400 300])
set(gca, 'fontsize', 20);
xlabel('Concentration (mM)')
ylabel('Max OD600')
box off;
xlim([0 125]);
